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Abstract 

We show that the explosion of the first supernovae can trigger low-mass star for- 
mation via gravitational fragmentation of the supernova-driven gas shell. If the 
shell mass does not exceed the host galaxy gas mass, all explosions with energies 
EsN ^ 10^^ erg can lead to shell fragmentation. However, the minimum ambient 
density required to induce such fragmentation is much larger, no > 300 cm~^, 
for Type II supernovae than for pair-instability ones, which can induce star forma- 
tion even at lower ambient densities. The typical mass of the unstable fragments 
is ~ IO^^^Mq; their density is in the range 110 — 6 x 10^ cm^'^. Fragments have 
a metallicity strictly lower than 10~^'^ Zq and large values of the gravitational-to- 
pressure force ratio / — 8. Based on these findings, we conclude that the second gen- 
eration of stars produced by such self-propagating star formation is predominantly 
constituted by low-mass, long-living, extremely metal-poor (or even metal-free, if 
mixing is suppressed) stars. We discuss the implications of such results for Pop III 
star formation scenarios and for the most iron-poor halo star HE0107-5240. 
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1 Introduction 

As many recent numerical (Bromm, Coppi & Larson 1999, 2002; Abel, Bryan 
& Norman 2000) and semi-analytical (Schneider et al. 2002, Omukai & In- 
utsuka 2002, Omukai & Palla 2003) studies have shown, the first luminous 
objects, the so called Population III (Pop III) stars, are likely to be very mas- 
sive. According to these studies, stars with a characteristic mass of 100-600 
Mq originate in the primordial gas, whereas low mass objects should not have 
formed due to the lack of heavy elements. High mass stars are indeed required 
to account for the unexplained Near Infrared Background (Salvaterra & Fer- 
rara 2003; Magliocchetti, Salvaterra & Ferrara 2003). Bromm et al. (2001) 
and Schneider et al. (2002) have shown that there exists a critical metallic- 
ity {Zcr — 10"^='=^ Zq) setting the transition from a high-mass to a low-mass 
fragmentation mode of star formation. The exact value of Zcr depends on the 
fraction of heavy elements that are depleted onto dust grains (Schneider et al. 
2003). In this scenario, a gas cloud with a metallicity Z ~ 10~^'^ Zq can lead 
to the formation of low mass stars if ~ 20% of the heavy elements is in dust 
grains (Schneider et al. 2003). The recent discovery of the most iron-poor star 
([Fe/H]= -5.3 ± 0.2) ever seen in our Galaxy (Christheb et al. 2002) could 
be an example of such low-mass metal poor stars originated from a gas with 
mean metallicity of 10^^'^ Zq (Schneider et al. 2003; but see Umeda & Nomoto 
(2003), who have pointed out that this star could have instead formed from a 
carbon-rich, iron-poor gas, with a corresponding mean metallicity of 10~^ Zq). 

An alternative potentially viable mechanism to form low-mass, metal-poor 
stars does exist. Indeed, Nakamura & Umemura (2001, 2002, hereafter NUOl 
and NU02) have studied the collapse and fragmentation of filamentary pri- 
mordial gas clouds using one- and two-dimensional hydrodynamical simula- 
tions coupled with the nonequilibrium processes of molecular hydrogen for- 
mation. They have shown that filaments with relatively low initial density 
{n < few X 10^ cm^'^) tend to fragment into dense clumps before the central 
density reaches 10^ — 10^ cm^^; the fragment mass is around 100 Mq. In con- 
trast, if a filament has initially a larger density, fragmentation continues until 
the clumps become optically thick to H2 lines, leading to a typical fragment 
mass of ~ 1 Mq. So the resulting IMF of first stars would be bimodal, with 
a low-mass peak around '^1 — 2 Mq and a high-mass peak at a few hundred 
Mq. Mackey et al. (2003) have suggested that these low-mass peak stars could 
form in the gas shocked by the explosion of the first generation of very massive 
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supernovae (SNe). 



In this paper, we put this hypothesis on a more quantitative basis. We will 
show that the shell is gravitationally unstable only for large explosion energies 
and that gravitationally unstable clumps are likely to form long-living, low- 
mass, metal-poor stars, that could be detected through dedicated surveys as 
the Hamburg/ESO objective prism survey (Christlieb et al. 2001). 

The paper is organized as follows. In Section 2 we calculate the SN-driven shell 
evolution and derive the conditions for gravitational instability. In Section 3 we 
examine the instability in the expanding shell for different explosion energies 
and interstellar medium (ISM) densities, whereas in Section 4 we calculated 
the properties of unstable fragments. In Section 5 we constrain the fragment 
metaUicity; Section 6 discusses the consequences for Pop III star formation 
scenarios. 



2 Evolution of the shell 



One can show that as far as the ambient gas pressure and cooling can be 
neglected, the shell evolution is described by the analytic expression derived 
by Sedov (1959): 



Rsnoci^Y't^f^ (1) 



where t is the time elapsed from the explosion, Esn is the total explosion 
energy, and po = fJ'TriHno is the density of the ISM (uq is the gas number 
density and n = 0.59 is the mean molecular weight of a primordial ionized 
H/He mixture). 

The further evolution of a SN-driven shell in the interstellar medium can be 
studied using the thin shell approximation (Ostriker & McKee 1988; Madau, 
Ferrara & Rees 2001). The momentum and energy conservation yield the fol- 
lowing relevant equations: 

j^iVshPoRsh) = 47ri?,\(n - P), (2) 



dt 



= -AnRi^PbRsh - VshniiMTb) , (3) 



where the subscripts 's/i' and '6' indicate shell and bubble quantities, respec- 
tively. Here, Rgh is the shell radius, Vgh = {ATT/3)Rlf^ is the volume enclosed by 
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the shell; the overdots represent time derivatives. P is the ISM pressure. The 
shell expansion is driven by the internal energy Ej, of the hot bubble gas, whose 
pressure is Pf, = Eh/27iRlj^ (for a gas with adiabatic index 7 = 5/3). Finally, 
n|^;,A(ff,) is the cooling rate per unit volume of the hot bubble gas, whose 
average hydrogen density and temperature are nH,b and 3],, respectively. The 
right-hand side of eq. (2) represents the momentum gained by the shell from 
the SN-shocked wind, while the right-hand side of cq. (3) describes the me- 
chanical energy input, the work done against the shell, and the energy losses 
due to radiation. 



As cooling becomes important for the swept up gas, the shock is no longer 
driven by the heated gas and the evolution enters the momentum conserving 
phase, satisfying the familiar 'Oort snowplow' solution, Rgh oc t^^^. 



2.1 Cooling time 



The cooling time is given by 



2 npsA.(Tpg^ 



where Upg = 720(7 + 1) /(7 ~ 1) is the post-shock density and Tps the post-shock 
temperature for an adiabatic shock, and A(Tps) is the cooling function for a 
primordial gas. The evolution of the bubble gas temperature is given by 
(Madau, Ferrara, Rees 2001) 

^ = 3^Rg, + ^PsH- (5) 

dt Rsh Psh 10 C2 RghPsh 



where Ci = lQ7rfimpri/25k and C2 = (1 25/39) vTyump and 77 = 6 x 10^'^ (cg.s. 
units) is the classical Spitzer thermal conduction coefficient (we have assumed 
a Coulomb logarithm equal to 30) . This relation closes the system of equations 
(2)-(3). 

If the shell lives for sufficiently long time (i.e. t,,„o/ <^ t), the gas cools down to 
a final temperature Tf ~ 300 K. This temperature is the minimum temper- 
ature provided by cooling from H2 molecules forming under nonequilibrium 
conditions in the post-shock gas (Ferrara 1998). 
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2.2 Perturbations in the shell 



The problem of the fragmentation of expanding shell has been adressed by 
several authors (Ostriker & Cowie 1981; Elmegreen 1994; Wiinsch & Palous 
2001; Ehlerova & Palous 2002). Here we apply these studies to the case of first 
SNe exploding in a medium of primordial composition. 



The linear growth of perturbations in an expanding shell was derived by 
Elmegreen (1994). Expansion tends to suppress (stabilize) the growth of den- 
sity perturbations owing to stretching, hence counteracting the self-gravity 
pull. The instantaneous maximum growth rate is 



3R 



sh 



+ 



2 / \ 2 

Rsh\ _^ I T^GpoRsh\ 



R 



■sh , 



3c 



s,sh 



1/2 



(6) 



where Cg^sk = {hTf/jimH)^^^ is the sound speed in the shell. Instability occurs 
only if a; > 0, or 



Rsh ^ 1 T^GpoRsh _ t 



J- ..^^u-'-'sn ■across /y\ 



where t cross ~ Rsh/cs,sh is the crossing time in the shell and t// is the free-fall 
time. Furthermore, if the fragment mass is close to (but slightly above) the 
Jeans mass, the a; > condition translates into Rsh/ Rsh < ^/tff- So, large 
shell velocity-to-radius ratios inhibit the formation of gravitationally unstable 
fragments. These relations are derived under the assumption of supersonic 
motion, i.e. Rghit) > Cs^, where Cgfl is the sound speed in the ISM. If the 
shock decays into a pressure wave before the onset of the instability (i.e. 
before a; > 0), no fragmentation will take place and the shell will be dispersed 
by random gas motions. 



A/ = 7T^[l + (l-err\ (8) 



The most rapidly growing mode has wavelength 
GpoR. 

where the dimensionless parameter ^ is given by 

'8'/^Rsh\ f7TGpoRsh\ 



Rsh J \ 3c, c/i 
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We will show later that the typical fragment size {Rf = A//2) is greater than 
the thickness of the shell. In the radiative phase the thickness is determined by 
the so-called cooling length, i.e. the distance travelled by the gas accelerated 
at the postshock velocity in a coohng time: ARgh ~ Rshtcooi- In this case it is 
likely that disk-like fragments form with mass 

Mf = npshARshR}, (9) 



where psh is the mean density in the thin shell, given by the ratio between its 
mass (equal to the shell swept up mass) and its volume Vsh 

4 Rsh 



For subsequent purposes we define the ratio of gravitational-to-pressure force 
evaluated in the direction along the fragment radius, i?/: 

^ ~ {3/2)kTf ■ ^ 



3 Instability in the expanding shell 



We solve numerically Eq. (2)-(3) and Eq. (5) for different values of the unper- 
turbed ISM density no, and of the explosion energy, Esn, following the time 
evolution of the shell radius (Rsh) and velocity (Rgh)- 

The unperturbed medium is assumed to be metal free and homogeneous. 
The SN progenitor ionizes the medium within its Stromgren radius, Rg ~ 
(3(5i?/47rnoQ;^)-'^/^, where Qh is the mass-dependent stellar ionizing photon 
rate and a is the hydrogen recombination rate to levels > 2. The radius of the 
ionized bubble reaches Rg is a recombination time tree = (ct^o)^^! however, 
the heated gas will start to freely expand at roughly its sound speed {cg w 10 
km/s for an ionized gas) if the dynamical time Rg/cg < i*, the lifetime of the 
star. In this case, by the time of the SN explosion the ionized region has grown 
to several tens of pc essentially independent of the initial value of Rg] inside 
this region the density and pressure are almost constant, thus motivating our 
constant density approximation. 

After the SN explosion, in the absence of ionizing photons, the gas would 

recombine on a timescale, tree- In practice, the gas is kept ionized by the radi- 
ation from the shock radiative precursor and its temperature remains ~ 10"^ K 
until ionizing photons continue to be produced. This condition is fuUfiUed as 
long as the shock velocity is larger than 80 km/s (Draine & McKee 1993). It is 
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only a recombination time after such event that the gas can finally recombine 
and cool to 300 K. 

At each time step we check the supersonic motion {Rsh > Cs,o) and instability 
{cu > 0) conditions. If both conditions are satisfied we further check that the 
cooling time is shorter than the age of the shell, so that the temperature has 
decreased to Tf ~ 300 K. At this point we calculate the fragment density (Eq. 
10), mass (Eq. 9), and the ratio / (Eq. 11). 

We consider here SN energies in the range 10^^ — 10^^ erg, the lower limit 
being typical of core-collapse Type 11 SNc (SNII) explosions and the upper 
limit corresponding to the energy released by the most massive pair- instability 
SN (SN^^) explosions (Heger & Woosley 2002). We explore the range of ISM 
densities no — 10~^ — 10^ cm~^. Fig. 1 shows the time-evolution of the shell 
radius and velocity for uq = 5 cm~^ and three explosion energies. The bottom 
panel shows the corresponding evolution of the instability growth rate cu. Only 
the explosions of SN^^ induce shell instability and fragmentation. 

Fig. 2 shows the region of the parameter space Esn — t^q where both the condi- 
tions a; > and Rsh > o are verified and therefore the shell is gravitationally 
unstable. In addition, we also require that the shell mass at the time of insta- 
bility does not exceed the total gas mass of the host galaxy. To exemplify, in 
the Figure we assume that the first stars form in halos with mass Mgai corre- 
sponding to virial temperature T^j^ — 10^ K. Under these conditions, we see 
from Fig. 2 that, if the shell mass docs not exceed the baryonic in the galaxy, 
all explosions with energies Esn > 10^^ erg can lead to shell fragmentation. 
However, the minimum ambient density required to induce such fragmenta- 
tion is much larger, uq > 300 cm~^, for SNII (explosion energy 10^^ erg) 
than for SN^^, for which the instability can occur even at lower ambient den- 
sities. Hence, pair-instability SNe are more suitable triggers of induced star 
formation. 

We conclude that the first pair- instability SNe are able to trigger self- propagating 
star formation under a wide range of ambient conditions, whereas expanding 
shells created by less energetic 'classical' or subluminous (as those proposed 
by Shigeyama & Tsujimoto 1998) SNII cannot fragment unless the density is 
very high. 



4 Properties of the fragments 

Let us now explore in more details the properties of the gravitationally unsta- 
ble fragments. As it is clear from Eqs. (8)-(9) the fragment mass is a function 
of the explosion energy (through Rg^ and Rsh), of the initial density of the 
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Fig. 1. Top panel: shell radius (increasing curves) and velocity as a function of time 

for three SN energies {solid line: Esn = 10^^ erg; dotted line: Esn = 10^^ erg; 
dashed line: Esn = 10^^ erg). The ISM density is no = 5 cm~^. Bottom panel: 
instability growth rate evolution for the three SN energies. The curves are plotted 
until the shock decays into a pressure wave. 

medium, uq, and of the density of the thin shell, Psh (proportional to the shell 
swept up mass, and thus function of Rgh)- The growth time of the unstable 
fragments is ~ o;"^. 

In all unstable cases, the age of the shell is larger than the cooling time and 
the thin shell approximation is valid {ARgh/ Rsh < 10"^). The instability sets 
in at 0.2-50 Myr after the explosion and the typical fragment radius and mass 
is in the range 2 — 1000 pc and 10'^"^ Mq, respectively (the upper limits refer 
to the lowest ambient density). The density of the fragments varies between 
110 cm~'^ and 6 x 10^ cm~^; as we discuss below, this spread has important 
implications for the mass of the second generation of stars formed. Large values 
of the gravitational-to-pressure force ratio (/ ~ 8) are found in all unstable 
cases. This is because high shell velocities stretch the perturbations, stabilizing 
the shell (see Eq. 7) and so the instability conditions are satisfied only when 
the fragments have collected a sufficiently large mass, yielding large values of 

/■ 
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Fig. 2. Region of the parameter space {Esn ~ ^t-o) where fragmentation occm's. See 
text for details. 

NUOl (see their Fig. 6) have shown that primordial filaments fragment into 
dense clumps whose masses depend on and /. For a quasi-equilibrium 
clump (i.e. / > 1), fragmentation proceeds until the fragment is close to the 
Bonnor-Ebert mass value corresponding to the initial density of the clump ^ 
(Palla 2002, Mackey et al. 2003). If / is increased to 3, a bifurcation takes place 
at ~ 10^ cm~^. For models with Ush ~ 10^ cm~^ the minimum fragment 
mass is 1-2 Mq, while for n^h ^ 10^ cm~^ it is larger than ~ 100 Mq. Further 
increase of / causes the low mass regime to extend down to ~ 10^ cm~^. 

In the low mass regime, the fragmentation proceeds down to smaller subclumps 
searching for an equilibrium between gravitational and pressure forces (thus 
decreasing /); at the same time the density of the subclumps tends to increase. 
For / values around 8, this equilibrium cannot be reached before the subclump 
becomes optically thick to H2 lines. At this stage, the initial clump has already 
fragmented into protostellar cores with typical masses ~ 10~^ ^ . Although 



^ This is strictly true only for n^h ^ 'f^cr = lO'^ cm~^ where ricr is the critical 
density which marks the transition from NLTE to LTE regime for H2 cooling (see 
discussion in Schneider et al. 2002). 

^ In the case of a strictly metal-free gas the minimum mass essentially coincides 
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the exact value of the resulting stellar mass depends on the details of subse- 
quent gas accretion and. possibly, on the merging of protostellar cores, it is 
unlikely that these are so efficient to form a massive star. In fact, the relax- 
ation time corresponding to the ensemble of m ~ 10~^ Mq protostellar cores 
resulting from the fragmentation of a Mf — 10^ Mq clump is ~ 36 Gyr. 

As mentioned above, wc find / 8 for all the unstable case. According to 
NU02 (see their Fig. 5a), this / value leads to the preferential formation of low- 
mass stars provided the density of the fragment is larger than ~ 2 x 10^ cm~^. 
Such condition is verified for all values of the ambient density Uq and explosion 
energy Esn except for the small area depicted in Fig. 2, where instead high 
mass (M 100 Mq) stars can form as a result of the fragment collapse. The 
formation of high mass stars can only be induced by SN^^, but the range of 
suitable ambient densities is small enough that such event can be regarded as 
unlikely. 



5 Mixing efficiency and metafiicity of the fragments 

In order to derive the mean metallicity of the thin shell (and of the fragments), 
one should be able to follow the mixing of SN ejecta with the shell gas. The 
standard hydrodynamical response to a sudden release of energy implies that 
the two gases (the swept up matter and the ejecta), after crossing their respec- 
tive shocks find themselves well separated by a contact discontinuity. However, 
many mechanisms (e.g. cloud crushing, thermal evaporation, hydrodynamical 
instabilities, effects caused by explosion inside wind-driven shells and by frag- 
mented ejecta) can disrupt the contact discontinuity leading to a mixing of 
the heavy elements into the swept-up matter in the thin shell (Tenorio-Tagle 
1996). A precise description of such physical processes is extremely difficult 
as it requires ultra-high resolution simulations (de Avillez & Mac Low 2002; 
Kifonidis et al. 2003). 

The mean metallicity of the shell (or fragment) is 

Zf = frrn.^^ (12) 

where fmix is the (unknown) mixing factor of ejected matter in the thin shell 
and Mej is the ejected mass in heavy elements. For SN-y-y, M^j — 0.4M* (Heger 

6 Woosley 2002). Even requiring that mixing is very efficient (i.e. fmix = 1), 
we always find that the mean metallicity of the fragments is 10~^-^ Zq < Z < 
10~^"^ Zq. In more realistic cases, however, mixing is likely to be much more 

with the Chandrasekhar mass, i.e. ~ 1 — 2Mq (Uehara et al. 1996) 
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inefficient {fmix <C 1), so the above value, 10~^'^ Zq is a strict upper limit. At 
face value, if the most iron poor star observed in our Galaxy ([Fe/H] = —5.3, 
Christlieb et al. 2002) formed through this mechanism, it would imply that 
fmix ~ 0.03. Note that in principle zero-metallicity low-mass stars can form 
via the proposed mechanism if fmi^ — 0. 



6 Discussion 

We have shown that the explosion of the first SNe are able to trigger further 
star formation. If the shell mass docs not exceed the baryonic in the galaxy, all 
explosions with energies Esn > 10^^ erg can lead to shell fragmentation. How- 
ever, the minimum ambient density required to induce such fragmentation is 
much larger, uq > 300 cm~^, for SNII (or subluminous supernovae) than for 
pair-instability ones, which can induce star formation even at lower ambient 
densities. Fragmenting shells become gravitationally unstable 0.2-50 Myr after 
the explosion, depending on the explosion energy and density of the unper- 
turbed medium. We have identified the fastest growing mode in the unstable 
shells and derived the mass, density and gravitational-to-pressure force ratio, 
/, of the corresponding fragments. Their typical mass is in the range 10^~^ Mq 
and density in the range 110 — 6 x 10^ cm~^. In all cases, we find very large 
values of the gravitational-to-pressure force ratio (/ ~ 8), indicating that the 
formation of low-mass, long-living stars is the most likely outcome. A narrow 
region of the Uq — Esn parameter space allows the formation of massive stars 
from fragments originating from SN^^ shells. However, the density range is 
small enough to consider this possibility as unlikely. 

The mean metallicity of the stars formed through this mechanism depends 
on the mixing history of the shell which is largely unknown, as discussed 
in the previous Section. Depending on the efficiency of this process and on 
the combined properties of the explosion energy and ambient density, the 
metallicity of this second generation, low-mass stars can be anywhere in the 
range < Z < lO'^-^ Zq. 

The mechanism here proposed is thus able to produce long-living, low-mass 
extremely metal-poor (or even metal-free, if fmix = 0) stars that can be found 
in the Milky Way halo. These stars can populate the low-mass peak of the 
bimodal IMF proposed by NUOl and NU02 or can be the typical members 
of the so-called II. 5 population (Mackey et al. 2003). The formation of these 
second generation stars requires the presence of a first generation of massive 
stars exploding as pair-instability SNe. 

So far, only one extremely iron-poor star with [Fe/H]< —5 has been identified 
(HE0107-5240, Christlieb et al. 2002). If this result were to be confirmed by the 
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analysis of the complete volume sampled by the Hamburg/ESO survey, several 
implications for the self-propagating star formation mode can be drawn: the 
observed lack of [Fe/H]< —5 halo stars would imply that mixing of metals in 
the unstable shell must occur with moderately high efficiencies {fmix <^ 10%). 
Alternatively, the mechanism here proposed might not lead to a large number 
of observed iron-poor halo stars because: (i) the typical mass of this second- 
generation stars is ^ 1 — 2Mq so that their lifetimes are not long enough to 
be observable as main sequence stars in our Galaxy halo; (ii) pair- instability 
SNe are rare because the IMF of the first stars is shaped so that only a 
small fraction of zero-metallicity massive stars form in the mass range of pair- 
instability SNe, UOMq ^ M^^ ^ 26OM0. 

Conversely, if more very iron-deficient ([Fe/H]< —5) stars were to be identified, 
the statistics and properties of these old stellar rehcs might lead to important 
constraints on the dominant processes which enable low-mass star formation 
in primordial environments. In particular, three viable mechanisms have been 
proposed for the origin of HE0107-5240 (Christlicb et al. 2002). Their main 
difference relies in the interpretation of the observed surface abundance of 
C, N (and O; Christheb et al. 2004). In particular, if one assumes that the 
observed Fe is not a good indicator of the metallicity of the gas cloud out 
of which the star formed and that C and N were already present in the star 
forming gas (C, N pre-formation scenario), then HEO 107-5240 formed with 
an initial metallicity of Z ~ 10~^Zq. Indeed, Umeda & Nomoto (2003) have 
shown that the abundance pattern of HE0107-5240 is in good agreement with 
nucleosynthetis yields of a faint SN explosion of a ~ 25Mq zero-metaUicity 
star releasing a kinetic energy of 0.3 x 10^^ erg. If so, due to the small ex- 
plosion energy, self-propagating star formation could not have occurred and 
the star must have formed on a longer timescale from the gas enriched by the 
SN. Alternatively, if one assumes that the observed Fe is a good indicator of 
the parent cloud metallicity and that C and N were synthetized in the stel- 
lar interior (C, N post-formation scenario) ^ , then HE0107-5240 formed with 
an initial metallicity of Z ^ Indeed, the observed abundance pat- 

tern for elements heavier than Mg (that cannot be formed in the interior of 
a O.8M0 star such as HE0107-5240 and thus retain memory of the nucleosyn- 
thesis yields of the pre-enriching star) are well reproduced by the predicted 
yields of a 200 M© SN^^. In this case, HE0107-5240 could have formed on 
a short timescale through the self-propagating star formation mechanism or, 
on a longer timescale, from the gas enriched by the SN if 20% of the metals 
were depleted onto dust grains (Schneider et al. 2003; Schneider, Ferrara & 
Salvaterra 2003). 



^ These "post-formation" scenarios require a sensible explanation for the large ob- 
served C and N abundances. So far, it is still unclear whether efficient mixing during 
the He-flash can fully account for the observed abundances or if it had enough time 
to operate in HE0107-5240 (Schlattl et al. 2001, 2002; Siess et al. 2002). 
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